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Abstract 

q ■ The Kuramoto model, which describes synchronization phenomena, is a system of 

CN ■ ordinary differential equations on Af-torus defined as coupled harmonic oscillators. The 

order parameter is often used to measure the degree of synchronization. In this paper, a 
few properties of the continuous model for the Kuramoto model are investigated. In par- 
ticular, the moments systems are introduced for both of the Kuramoto model and its con- 
CN ■ tinuous model. By using them, it is proved that the order parameter of the iV-dimensional 

Kuramoto model converges to that of the continuous model as N — » oo. 

Q 

43 : 1 Introduction 

Collective synchronization phenomena are observed in a variety of areas such as chemical 
reactions, engineering circuits and biological populations lfT8l. In order to investigate such 
a phenomenon, Kuramoto lfTTl[T2l proposed the system of ordinary differential equations 

(N . de- k ^ 

O ■ — i = cot + - V sin(0,- - Bd, i = 1, • • • ,N, (1.1) 

IT) 1 at N 

7=1 

m 

where t e [0, 2n) denotes the phase of an z'-fh oscillator on a circle, a», 6 R denotes 
its natural frequency, K > is the coupling strength, and where Af is the number of 
oscillators. Eq. dl.ll) is derived by means of the averaging method from coupled dynamical 
systems having limit cycles, and now it is called the Kuramoto model. 

It is obvious that when K = 0, 6i(t) and 8j(t) rotate on a circle at different velocities 
unless u>i is equal to ojj, and it is true for sufficiently small K > 0. On the other hand, if K 
is sufficiently large, it is numerically observed that some of oscillators or all of them tend 
to rotate at the same velocity on average, which is called the synchronization [HH I2TTI . 
If N is small, such a transition from de-synchronization to synchronization may be well 
revealed by means of the bifurcation theory |[6l \T3[ O. However, if is large, it is 
difficult to investigate the transition from the view point of the bifurcation theory and it is 
still far from understood. 
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In order to evaluate whether synchronization occurs or not, Kuramoto introduced the 



which gives the centroid of oscillators, where r, iff e R. It seems that if synchronous state 
is formed, r(f) takes a positive number, while if de- synchronization is stable, r(t) is zero 
on time average. 

The infinite-dimensional version or the continuous version of the Kuramoto model 
have been well investigated to reveal a bifurcation diagram of the order parameter (see 
ttH El SI d [13 [IS EU E21 E3B and references therein). Such systems are rather tractable 
because the order parameter for the infinite-dimensional version can be constant in time, 
while the order parameter for the finite dimensional Kuramoto model is not constant in 
general because solutions fluctuate due to effects of finiteness ETI . Now the questions 
arise : How close is the order parameter of the infinite-dimensional version to that of the 
finite-dimensional Kuramoto model? What is the influence of finite size effects? This 
issue has been studied by many authors, see a reference paper [TJ by Acebron et al. In 
particular, Daido ^ found the scaling law \r-(r N (t))\ ~ (K c -K)~ l/2 N~ l/2 forK < K c and 
\r - (r N (t))\ ~ (K — K c )~ l/s N~ l/2 for K > K c , although his analysis is not rigorous from a 
view point of mathematics, where r is the order parameter of the infinite-dimensional ver- 
sion in a steady state, r N (t) is the order parameter of the iV-dimensional Kuramoto model, 
( ) denotes the time average, and where K c is the transition point from de-synchronization 
to synchronization. 

In this paper, the continuous model, the continuous limit of the Kuramoto model, is 
introduced. It is proved that the order parameter of the iV-dimensional Kuramoto model 
converges to that of the continuous model in the sense of probability, and their difference is 
of 0(N~ 1/2 ) as N — » oo for each t (note that we do not take time average). To prove this, the 
(m, fc)-th moments are defined for both of the continuous model and the Kuramoto model. 
In particular, (0, l)-th moment is the Kuramoto's order parameter. It is remarkable that 
both of the continuous model and the Af-dimensional Kuramoto model become the same 
evolution equation, called the moments system, if they are rewritten by using the moments. 
It means that any solutions of the continuous model and the Af-dimensional Kuramoto 
model for any N are embedded in the phase space of the moments system. This fact 
allows us to measure the distance between a solution of the continuous model and that of 
the iV-dimensional Kuramoto model in the same phase space. These results and the central 
limit theorem prove that the difference between the order parameter of the iV-dimensional 
Kuramoto model and that of the continuous model is of order 0(N~ 1/2 ) for each t, provided 
that initial values and natural frequencies for the A^-dimensional Kuramoto model are 
independent and identically distributed according to a suitable probability measure. 



order parameter r{t)e by 




(1.2) 
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More generally, a globally coupled phase oscillators defined to be 
dQ K N 

i =Wi + ^Z /( ^'" 9i) ' i=h---,N, (1.3) 

is called the Kuramoto-Daido model D9JI21, where the 2^-periodic function / : R — > R is 
called the coupling function. The results in this paper are easily extended to the Kuramoto- 
Daido model. 



2 Continuous model 



In this section, we introduce a continuous model of the Kuramoto model and show exis- 
tence, uniqueness and other properties of solutions of the model used in a later section. 

Let us consider the Kuramoto model (11.11) . Following Kuramoto, we introduce the 
order parameter Z° by 



(2.1) 



The quantities ZT will be defined in the next section. By using it, Eq. dl.ll) is rewritten as 



ddi 
dt 



K 



OJi + 



2V-T 



(Z\{t)e 



(2.2) 



where Z° denotes the complex conjugate of Z°. Motivated by these equations, we intro- 
duce a continuous model of the Kuramoto model, which is an evolution equation of a 
probability measure p, = p t (6, to) on S 1 x R parameterized by t e R, as 



Z?(0 := f f 
JrJo 

p (6, to) = h(9, to), 



e^°dp t , 



(2.3) 



where h{0, to) is an initial measure. The Z^(t) is a continuous version of Z^(f), and we also 
call it the order parameter. If we regard 

K 



v := to + 



as a velocity field, Eq. fl2.3l) provides an equation of continuity dpjdt + d(p t v)/d8 = 
known in fluid dynamics. It is easy to prove the low of conservation of mass: 



\ \ XE(co)dp t = I I XE(oj)dh =: g(E), 
JrJo JrJo 



(2.4) 
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where E is any Borel set on R and / f £ (6t») is the characteristic function on E. A function 
g defined as above gives a probability measure for natural frequencies co e R such that 
f R dg = 1. In particular ^dp t = 1 if fj^dh = 1. 

By using the characteristic curve method, Eq. (l2.3l) is formally integrated as follows: 
Consider the equation 

^ = w + -^=(Z?(f)c- V=rx -zJa)c V=Tx ), xeS 1 , (2.5) 
dt 2V-T 

which defines a characteristic curve. Let x = x(t, s; 9, co) be a solution of Eq. (l2.5l) satis- 
fying the initial condition x{s, s; 9, co) = 9. Then, along a characteristic curve, p, is given 
as 

p t (6, co) = h(x(0, t; 9, co), co) exp[^ f {7P x {s)e' ^~ U(sM + ^is)e^ x(s ' tM )ds\. (2.6) 

2 Jo 

By using Eq. (l2.6l) . it is easy to show the equality 

J J a{9, co)dp t = J J a(x(£, 0; 0, o>), co)dh, (2.7) 
JrJo JrJo 



for any measurable function a{9, co). In particular, the order parameter Zj '(f) are rewritten 

as 

e^^^dh. (2.8) 



z?(o = f r % 

JrJo 



Substituting it into Eqs. (l2.5l) . (12.61) . we obtain 

d r r 2n 

—x(t, s;9,co) = co + K J J sin(x(f, 0; 0', oj') - x(f, 5; 0, co))dh(9', co'), (2.9) 

and 

ds- I I cos^x(s,0;9',co')-x(s,t;9,cofjdh(9',co')j. 

(2.10) 

Even if h(9,co) is not differentiable, we consider Eq. (l2. 101) to be a weak solution of 
Eq. fl2.3l) . Indeed, even if h and p, are not differentiable, the quantity (12.71) is differen- 
tiable with respect to t when a(9, co) is differentiable. 

Theorem 2.1. (i) There exists a unique weak solution p t of the initial value problem 



231) for any f > 0. 



(ii) Solutions of (12.31) depend continuously on initial measures with respect to the weak 
topology in the sense that for any numbers T, s > and for any continuous function 
a(9, co) on S 1 x R, there exist numbers M(b) > and 5 = 8{T, s, a) > such that if initial 
measures h\,h 2 satisfy 

b{9, co)(dh[ - dh-i) 



If 

JrJo 



<M(b)S, (2.11) 
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for any continuous function b(8, to), then solutions p t l and p fj2 with p 0l =h\ and p ,2 = ^2 
satisfy 



r r ace, 

JrJo 



toXdpt,i - dp u2 ) 



< e. 



(2.12) 



for < ? < r. In particular if a is Lipschitz continuous, then s ~ 0(8) as 6 — » 0. 

Proof of (i). It is sufficient to prove that the integro-ODE (12.91) has a unique solution 
x(t, s; 6, to) satisfying x(s, s; 8, to) = 8 for any t, s > and 8 e S 1 . Let us define a sequence 
{x n (f,O;0,w)}~ o tobe 



-2n 



x n+ i(t,0;8,a>) = xo(t,0;8,a>) + K I dr - I J f(x n (r, 0; 0', a/) - x„(r, 0; 8, a>))dh(8', to') 

Jo JrJo 

(2.13) 

and xq(?, 0; 8, to) = 8 + o>f, where /(0) = sin 8 (since we prove the theorem for any C 1 
function f(8), the theorem is also true for the continuous model for the Kuramoto-Daido 



model (11.31) ). We estimate \x n +\(t, 0; 8, to) - x n (t, 0; 8, to)\ as 
\x„+i(t, 0; 8, to) - x n (t, 0; 8, to)\ 

■*2n 



< K I dr ■ 

Jo 



If 

JrJo 



< 



f(x n (T,0;d / ,to')-x n (r,0;d,to)) 

-f(x n -\(j, 0; &, to') - x„_i(t, 0; 8, to)) dh(ff , to') 

KLldT-l (|x„(t,O;^V)-x„-i(t,O;0W)| 
Jo JrJo 

+|x«(t, 0; 8, to) - x„-i(t, 0; 8, to)\}dh(6', to'), 



where L > is the Lipschitz constant of the function f(8). When n = 0, we obtain 



|jci(f,O;0,aO-;co(f,O;0,aO| < K 

< KMt, 



Jo JrJo 



|/(x (t, 0; 0', w') - x (r, 0; 8, to))\ dh(8', to') 



where M = max \f(8)\. Thus we can show by induction that 

\x n (t,0;8,to)-x n . l (t,0;8,to)\ < 2 n ~ l K n M— . (2.14) 

n\ 

This proves that x n (t, 0; 8, to) converges to a solution of the equation 
dx C C 2n 

—(t,0;8,(o) = a> + K f(x(t,0;8',to') - x(t,Q;8,to))dh(8' ,to'), 

at J R Jo 

as n — » 00 for small f > 0. Existence of global solutions are easily obtained by a standard 
way because the phase space S 1 is compact; that is, solutions are extended for any t > 0. 
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Uniqueness of solutions is also proved in a standard way and the detail is omitted. With 
this x{t, 0; 9, to), we define a sequence {x n {t, s; 9, to)}™ =0 to be 



x n+1 (t, s; 9, to) = x (t, s; 9,to) + K 



Jo JrJo 



f{x(r, 0; ' , to') - x n (r, s; 9, to))dh(9 , to ) 

(2.15) 

and xo(t,s;9, to) = 9 + to{t - s). Then, existence and uniqueness of global solutions 
x(t, s; 9, to) is proved in the same way as above. For this x(t, s; 9, to), Eq. (l2. 101) gives a 
(weak) solution of Eq. (l2.3l) . 

Proof of (ii). Suppose that initial measures h\, h 2 satisfy Eq. (l2.1 II) . Let p t ^ and p u2 be 
solutions of Eq. fl2.3l) satisfying p ,i = h\ and p ,2 = h 2 . Let x { = x t (t, 0; 9, to), (i = 1, 2) be 
solutions of 

f f f(Xi(t, 0; 9', to') - x,(t, 0; 9, to))dhi(9', to'), x,eS l , (2.16) 
JrJo 



d -^ = co + K 

dt 



satisfying x,(0, 0; 9, to) = 9, respectively. Then we obtain 



— (*ia,0;6>, to)-x 2 (t,0;9, to)) 
dt 



= K 



= K 



f f / '(jcia,0; 9', to') -X!(t,0; 9,oo))dh x {& , to') 
JrJo 

-K f f f(x 2 (t,0;9',to')-x 2 (t,0;9,to))dh 2 (9',to') 
JrJo 

f f /(jCi^O^'^O-Xi^O;^^))^!^',^')-^^',^')) 
JrJo 

+K I Io (^ l(r ' 0; ^' w ' ) - Xl(? ' 0; ^ £t>)) 

-/(x 2 (f, 0; 0', to') - x 2 (t, 0; 0, to)))dh 2 (9', to'). 

Integrating it yields 

\xi(t, 0; 9, to) - x 2 (t, 0; 9, to)\ 

< \ K \ \ /(xi^O^'^O-^i^O;^^))^!^',^')-^^',^')) 
Jo JrJo 



(2.17) 



ds 



+ 



Jo JrJo 



f(xi(s, 0; & ', to') - xi(s, 0; 9, to)) 

dh 2 {9',to')ds 



-f(x 2 (s,0;9',to') - x 2 (s,0;9,to)) 

(\xi(s,0;ff ',to') - x 2 (s,0;9' ',to')\ 
+\xi(s, 0; 9, to) - x 2 (s, 0; 9, to)\}dh 2 (9', to')ds, 



(2.18) 



where L is the Lipschitz constant of / and M is a constant arising from Eq. (l2.1 II) . If we 
put 



F(t, 9, co) = |*i (f, 0; 9, co) - x 2 (t, 0; 9, co)\, 



then (12.181) provides 

f f F(t,9,co)dh 2 <KMSt + 2KL 
JrJo 

Now the Gronwall inequality proves 

r c ln „, „ „ Mi 

JrJo 

Substituting it into (12.181) yields 



r r 

Jo JrJo 



F(s, 9, to)dh 2 ds. 



F(t,9,oj)dh 2 <'-^-(e 2KLt -l). 



|*i(f,O;0, u))-x 2 (t,0;9, co)\ 



< KMSt + 



KM6 



\)ds + KL \x\(s, 0; 8, co) - x 2 (s, 0; 8, co)\ds 



KMSt M8 
< : — + — (e 



f (e 2KLs - l)ds + KL f 
Jo Jo 

f 

Jo 



2 4L 

The Gronwall inequality is applied again to obtain 



1) + KL \xi(s,0;6,(o)-x 2 (s,0;6,(o)\ds. 



M6 



\xi(t, 0; 9, co) - x 2 (t, 0; 9, co)\ < '-^(e 2KLt - D- 



(2.19) 



Finally, the left hand side of Eq. (l2.12l) is estimated as 

I I a(9, to)(dp trl - dp ta ) 
JrJo 

= J" ^ {a(x\(t,Q;9,co),co)dh\ - a(x 2 (t,0;9,co),co)dh 2 j 

< I I \a{x\{t, 0; 9, co), co) - a(x 2 (t, 0; 9, co),co)\ dh 2 
JrJo 

+ J I a(xi(t,0;9,co),co)(dhx - dh 2 
JrJo 



(2.20) 



Since a(9, co) is continuous and since Eq. (12. 191) holds, the first term in the right hand side 
of the above is less than e/2 for < t < T if 5 is sufficiently small. The second term is 
also less than e/2 if 8 is sufficiently small because of Eq. (l2.11l) . This proves Eq. (l2.12l) . It 
is easy to see by Eq. fl2.20l) that if a(9, co) is Lipschitz continuous, then s is of order 0(6). 
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3 Moments system 



In this section, we introduce a moments system to transform the finite-dimensional Ku- 
ramoto model (11.11) and its continuous model (12.31) into the same system. We prove by 



using the moments system that the order parameter (11.21) for the Kuramoto model con- 
verges to the order parameter Z°(Y) for the continuous model as N — » oo under appropriate 
assumptions. 

For a given probability measure h{6, w) onS'xR, suppose that absolute moments 

M n k := f f We^^ldh (3.1) 
JrJo 

exist for k = 0, +1, • • • and m = 0, 1, • • • . Then, the moments m' k are defined to be 

ml ■= f f aj n e^ ko dh. (3.2) 
JrJo 

Conversely, if there exists a unique probability measure h for a given sequence of numbers 
i m " k }n,k such that Eq. (l3.2l) holds, then h is called M-determinate. In this case, we also say 
that moments {m n k } n ^ is M-determinate. Many (necessary and) sufficient conditions for 
which h is M-determinate have been well studied as the moment problem [|2] [191 OH 
l20ll . For example, one of the most convenient conditions is that if h has all absolute 
moments M' k and they satisfy Z^iC^o + l)~ 1/ ' n = 00 (Carleman's condition), then h is 
M-determinate. 

In what follows, we suppose that an initial measure h(9, oS) for the initial value prob- 
lem (12.31) has all absolute moments and is M-determinate. 

Example 3.1. If h has compact support, then h is M-determinate. Suppose that h has 
a probability density function of the form h(6)g(co). If g(oj) is the Gaussian distribution, 
then h is M-determinate. If g(oj) = 1/(^(1 + co 2 )) is the Lorentzian distribution, h is not 
M-determinate because h does not have all moments. 

Recall that a probability measure g for the natural frequency a> is defined through 
Eq. fl2.4l) . Since h has absolute moments 



M k 



» = f f \afe^ k6 \dh = f \(o\ n dg < oo, (3.3) 
JrJo Jr 

g also has all moments //„ := J R a) n dg, n = 0, 1, 2, • • • . Consider the Lebesgue space 
L 2 (R,dg). Since all moments ix m of g exist, we can construct a complete orthonor- 
mal system {P m (co)}°° =() on L 2 (R, Jg), by using the Gram-Schmidt orthogonalization from 



>m=Q i 



>m=0 

such that 



(Pn, PJ = £Pn(.C0)P m (C0)dg = j J ^ ^ (3.4) 
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where ( , ) denotes the inner product on L 2 (R, dg) and P n (co) is a polynomial of degree n. 
In particular, Po(co) = 1. It is well known that P n (co) satisfies the relation 



coP n (to) = b (to) + a n P n (io) + b 

for n = 0, 1, 2 • • • , where a n and b n are real constants determined by g 
defined as 

f a b 



M 



bo a x by 
b\ a 2 b 2 



(3.5) 
. The matrix M 

(3.6) 



is called the Jacobi matrix for g. Eq. (l3.5l) shows that the Jacobi matrix gives the / 2 (Z> ) 
representation of the multiplication operator 



M : p(co) i — > cop(co) 
on L 2 (R, dg), where Z 2 (Z> ) = {{x n }~ | Z,T=o l^l 2 < °°}- 



(3.7) 



If an initial measure h is M-determinate, so is a solution p, of the continuous model 
(12.31) because of Eq. (12.10l) . Let us define the (m, k)-th moments for p, to be 

Z?(f)= [ [ P m (io)e^ ke dp t , (3.8) 
JrJo 

for m = 0, 1, 2 • • • and = 0, ±1, +2, • • • . In particular Z^(t) is the order parameter given 
in Eq.d273]), and ZTj® = Z'"(t). Note that 



1 (m = 0) 
Z()W " l (m^O) 



are constants. It is easy to verify that 



K(t)\ < i 



(3.9) 



(3.10) 



by using the Schwarz inequality. By using Eq. (l2.7l) . an evolution equation for Zf (t) is 
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derived as 



dt 



P m (co)e^ kxme ^dh 



= - rr 

dt J R Jo 
JrJo 

= v-r* r r uPMe^^^dh 

JrJo 
2 JrJo 

= ^P\k J (b m P m+ i(a)) + a m P m ((jS) + bm-iPn-iia^e^^dpt 

+^ f r J P m (o;)(Z?(r) e ^- 1 ^-Z 1 (r) e ^ +1 ^)rf A 
^ JrJo 



<^\k{b m Z^ x + a m Z™ + b^ZT 1 ) + -rlA^lx ~ ^ti) 



(3.11) 



Put Zk = (Z?,Zl,Z%, ■ ■ ■ ) T , where T denotes the transpose. Define the Jacobi matrix M 
and the projection matrix P to be Eq. (l3.6l) and 



P = 



' 1 




(3.12) 



respectively. Then, Eq. d3.lll) is rewritten as 



d_ 

dt 



Z 2 

z 3 

v : / 



^P\M + -P 

2 



2V-TM 

3 V^M 



z 2 
z 3 



K 
+ 2 



-z\ Z 2 

2(Z°Z 1 -Z° 1 Z 3 ) 
\Z\Z 2 -Z\Z A ) 



(3.13) 



Note that equations for Z_i,Z_ 2 , • • • are omitted because Z_* = Z^,. The first term is a 
linear term and the second is a nonlinear term. We call Eq. (l3. Ill) or Eq. (l3.13l) the moments 
system. The dynamics of the system will be investigated in forthcoming papers 010. 



Let MD be the set of M-determinate sequences {ZT} m k in the sense that if {Z! 



k >>n,k 



MD, then there exists a unique measure h on S 1 x R such that Z' k n = f R £' P m ((o)e ^~ lke dh. 
Since there is a one-to-one correspondence between elements of MD and M-determinate 
measures, Thm.2.1 is restated as follows. 
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Theorem 3.2. (i) There exists a unique solution {Z' k (t)} mk e MD of the moments system 
if an initial condition {Z' k (0)} m ^ is in MD. 

(ii) Let (Zf(f)} and {Zf(t)} be solutions of the moments system with initial conditions 
}Z™(0)}, {Z' k (0)} e MD, respectively. For any positive numbers T and s, there exist posi- 
tive numbers C m ^ and 6 = 6(T, s) such that if 

|Z™(0)-Z™(0)|<C m ,,5, (3.14) 

for any m, k, then the inequality 

\Z?(t)-Z?(t)\<e (3.15) 
holds for < t < T. In particular s ~ 0(6) as 6 — * 0. 



For the yV-dimensional Kuramoto model (11.11) . we define the (m, k)-th moments to be 



ZHO-^^K)^^, (3.16) 

7-1 

for m = 0, 1 , 2, • • • and £ = 0, ±1 , +2, • • • . In particular Z° is the order parameter defined in 
Eq. (l2.1l) . By using Eq. dl.ll) . it is easy to verify that Z n k (t)\ satisfy a system of differential 
equations 

d^k W J t, ( , A m+ i 7in. u 7 m-\\ , kK 707 - 

— = y-ik\b m Z k +a m Z k +b m _iZ k ) + —(Z l Z k _ l - Z^Z k+l ). (3.17) 

It is remarkable that Eq. (l3.17l) has the same form as Eq. (l3.1 II) . This means that all solu- 
tions of the Kuramoto model for any yV are embedded in the phase space of the moments 
system (13.1 II) . This fact allows us to prove Theorem 3.3 below. Originally the moments 
system for the Kuramoto model was introduced by Perez and Ritort IfTTII . although their 
definition of the moments is H' k n := 1/7Y ■ Yu N j=\ <jj'Je^ k9i(t) . Since we adopt orthogonal 
polynomials {P m (w)}^ =0 to define moments (13.161) . our moments system is more suitable 
for mathematical analysis JH[5]|. 

Now we are in a position to show the main theorems in this paper, which states that 
differences between moments Z' k (t) and Z k (t) are of 0(1/ VAO and thus the continuous 
model (12.31) is proper to investigate the Kuramoto model (11.11) for large N. 

Theorem 3.3. Let p t be a solution of the continuous model (12.31) such that an initial 
measure h(6, cS) is M-determinate. Suppose that for the yV-dimensional Kuramoto model 
(11.11) . pairs (0j(O), loj) of initial values 8j(0), j = 1, • • • ,N and natural frequencies ojj, j = 
1, • • • , N are independent and identically distributed according to the probability measure 
hifi, to). Then, moments Z*(f) and 2*(t) defined by Eqs.dXl and (I3T61) satisfy 

|Z"(0-Zf(0l^0, a.s. (N^oo), (3.18) 
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for any m, k and t (a.s. denotes " almost surely"). Further, for any positive number 6, there 
exists a number C = C(m, k, t, 6) > such that 

P( \Z'l\t) - Zf (01 < CI Vtf ) -» 1 - 5, (V -> oo), (3.19) 

where P(A) is the probability that an event A will occur. 
Proof. The average of Zf(0) is calculated as 

E[Zf(0)] = EE^P.^^^] 

;'=i 

= E[P m (a>j)e^ ke ^] 

In 

PJ^e^dh = Zf(0). (3.20) 

Thus Eqs. (l3. 1 81) and (13.191) for t = immediately follow from the strong law of large 
number and the central limit theorem, respectively. Note that the strong law of large num- 
ber and the central limit theorem are no longer applicable for t > because 9j(tYs are 
not independent and identically distributed when t is positive. However, since Zf (?) and 
Zf (t) satisfy the same moments system, and since solutions of the moments system are 
continuous with respect to initial values (Thm.3.2 (ii)), Eqs. (l3.18l) . (l3.19l) hold for each 
positive t if they are true for t = 0. ■ 
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